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Abstract 

The tunneling effect of two-dimensional Dirac fermions in a constant magnetic field is studied. 
This can be done by using the continuity equation at some points to determine the corresponding 
reflexion and transmission coefficients. For this, we consider a system made of graphene as super- 
position of two different regions where the second is characterized by an energy gap t'. In fact, 
we treat concrete systems to practically give two illustrations: barrier and diode. For each case, 
we discuss the transmission in terms of the ratio of the energy conservation and t'. Moreover, we 
analyze the resonant tunneling by introducing a scalar Lorentz potential where it is shown that a 
total transmission is possible. 
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1 Introduction 



After the experimental realization of graphene in 2005, this kind of system became an attractive 
subject not only for experimentalists but also for theoretician physicists. This is because of the nature 
of its structures and the behavior of its relativistic particles. In addition of exhibiting the anomalous 
quantum Hall effect [TJ [2] , graphene gives an example of condensed matter physics where the quantum 
electrodynamics tools can be applied [3]. These new developments offered a laboratory for many 
investigations where interesting results are obtained by solving different problems. Because of the 
relativistic nature of their fermions, the system made of graphene renewed the interest of studying the 
Dirac fermions in two-dimensions. 

One of characteristics of Dirac fermions in graphene is their ability to tunnel through a potential 
barrier with probability one [U [5] . This so called Klein tunneling of chiral particles has long ago been 
proposed in the framework of quantum electrodynamics [SIEIIH], but was never observed experimen- 
tally. As appealing as the Klein tunneling may sound from the point of view of fundamental research, 
its presence in graphene is unwanted when it comes to applications of graphene to nanoelectronics. 
This comes about because the pinch-off the field effect transistor may be very ineffective. The same 
may occur because of the minimum conductivity of graphene at the neutrality point. One way to 
overcome these difficulties is by generating a gap in the spectrum. From the point of view of Dirac 
fermions this is equivalent to the generation of a mass term. 

A possibility of generating gaps in the graphene spectrum is to deposit graphene on top of hexagonal 
boron nitride (BN) [9]. This material is a band gap insulator with a boron to nitrogen distance of the 
order of 1.45 A [T2] (in graphene the carbon-carbon distance is 1.42 A) and a gap of the order of 4 eV. 
It was shown that in the most stable configuration, where a carbon is on top of a boron and the other 
carbon in the unit cell is centered above a BN ring, the value of the induced gap is of the order of 53 
meV. Depositing graphene on a metal surface with a BN buffer layer leads to n— doped graphene with 
an energy gap of 0.5 eV [10J. 

Theoretically, the tunneling effect of system type Si02-BN in zero magnetic field is discussed [TT] . 
In fact, it is assumed that it is possible to manufacture slabs with Si02-BN interfaces, on top of 
which a graphene flake is deposit. This will induce spatial regions where graphene has a vanishing gap 
intercalated with regions where the BN will cause a finite gap. The graphene physics is considered 
in two different regions: the k— region, where the graphene sheet is standing on top of Si02, and a 
q— region, where a mass-like term is present, caused by BN, inducing an energy gap of value 2t' . The 
effect of chiral electrons in graphene through a region is studied where the electronic spectrum changes 
from the usual linear dispersion to a hyperbolic dispersion, due to the presence of a gap. It is shown 
that contrary to the tunneling through a potential barrier, the transmission of electrons is, in this 
case, smaller than one for normal incidence. 

Motivated by the reason discussed above and in particular the investigation made in [11] , we would 
like to reply an interesting question, the fact that what happens to the tunneling effect of Si02-BN in 
the presence of an external magnetic field B. Otherwise, still we have the same conclusions reached 
for the case B = in [11] or the will be affected. To answer these inquiries, we study such system in 
the same conditions as in |llj but taking into account the effect of the gauge field. This will allow us 
to deal with some issues and end up with different conclusions. Moreover, under some conditions we 



1 



show that there is possibility for a total transmission. 

More precisely, we consider a system composed of two different regions, where the second is char- 
acterized by an energy gap t', in a perpendicular magnetic field. The energy spectrum solutions are 
obtained for both regions in terms of two Landau levels and t' . From the energy conservation E, we 
obtain a set of the energy values that allows us to discuss the tunneling effect By inspecting these 
values, one can end up with three limiting cases, which they have interesting consequences on the 
reflexion and transmission of the present system. 

To be concrete, we give two examples of system made of graphene. In the beginning, we consider a 
barrier in magnetic field and study the tunneling effect. Indeed, from the continuity equation, we get 
different solutions, which allowed us to explicitly determine the reflexion and transmission coefficients 
for different regions. They are used to show that the probability condition is satisfied, namely the 
sum of these coefficients is one. To characterize the transmission behavior we give different figures, 
which underline its properties in terms of the energy ratio y. Using the three limiting cases of energy 
we discuss the possibility to obtain a total transmission and give different interpretations. 

As far as the diode in magnetic field is concerned, we discuss the tunneling effect by splitting 
the whole system in three regions where the second different from the first and third, which they are 
identical. After getting different coefficients, we use an appropriate definition to show that the sum of 
reflexion and transmission is one. Analyzing this under some conditions, we conclude that it is possible 
to obtain a total transmission. We also discuss a limiting situation where the barriers are described 
by a scalar Lorentz potential. Using the boundary conditions, we derive different coefficients, which 
lead again to verify the probability condition and give different discussions about such potential. 

The present paper is organized as follows. In section 2, we formulate our problem by writing 
the Hamiltonian's describing two regions of our system. The eigenvalue problems will be solved to 
obtain the energy spectrum and its eigenspinors. We inspect the eigenvalue solutions from the energy 
conservation point of view to underline their properties in section 3. After establishing all needed 
materials, we treat the first illustration of our system, i.e. the barrier, in section 4. In section 5, 
we consider the diode in magnetic field as a second illustration where an interesting limit will be 
investigated. We conclude and give some perspectives in last section. 

2 Hamiltonian formalisms 

As claimed before, our system is a superposition of two different regions (I) and (II), with (II) has an 
energy gap t' . For this, we start by setting the necessary tools needed to treat each region separately. 
These concerns to write the corresponding Hamiltonian's and determine their eigenvalue solution as 
well as the eigenspinors. 

Region (I) can be identified to a two-dimensional subsystem of Dirac fermions where the Hamil- 
tonian for one massless relativistic fermion in the presence of a perpendicular magnetic field is 



VF <T ■ 7T 



(1) 



where vf ~ 10 6 



ms 1 is the Fermi velocity and a = (a x , a y ) are the Pauli matrices 




(2) 
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The conjugate momentum is n = p + ~ A, with A is a gauge field. Recall that, region (II) is of finite 
mass due to tf, then its appropriate Hamiltonian can be written as 



= VF d.jf + t'a z . 



(3) 



Its clear that, the mass term t'a z makes difference between both the above Hamiltonian's. This will 
play a crucial role in the forthcoming analysis. 



2.1 First region spectrum 

To do our task, we start by determining the energy spectrum and its eigenspinors for each involved 
region. In doing so, let us start by writing the Hamiltonian ([TJ as 



v F 





ir x + m. 



(4) 



Choosing the Landau gauge A = B (0, x,0), the momenta components take the form tt x = p x and 
% = Py + ^c~ x - These can be used to maps into 



( Px-iPy- i^r-x 

VF R 

\ Px + ip y + i^-x 







(5) 



As usual to get the energy solutions of © , one can use the eigenvalue equation for a given spinor 
0« = I ^ J of if®. This is 



\ <P2 



It implies the "kinetic balance" relation 



V f2 



eB 



"V ( '/',• +Py + — •'• ] ^2 =■ t ' ■ . 



eB 



IV F [ ~ip x +Py + —X ] (fl = ^'lf 2 - 



(6) 

(7) 
(8) 



To determine these spinor components, we can map for instance ([7|) into ([8]) to obtain a Schrodinger- 
like equation for ip2- After calculation, we obtain 
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> / eB 
Px + 1 Py + ~ x 



eBH 



<P2 



(9) 



It is similar to that of the harmonic oscillator up to some constant term and thus ip2 can be seen as its 
eigenfunction. To clarify this statement, let us consider the Landau Hamiltonian in the same gauge, 
such as 



1 

2m 



eB 

Px + I Py H x 



One can easily show that its energy spectrum reads as 

E®' = hw c (n + , n = 0,1,2, ■ 



(10) 



(11) 
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where the cyclotron frequency u> = In terms of the cylindrical parabolic functions, The corre- 
sponding wavefunctions are given by 



(x,y) = D n (Q)e lk yy 



(12) 



where k y is a wave vector along y-direction in region (I). We have set Q = x ~[^° and xq = k y l B , with 



the magnetic length Ib = \/ jr- It is not hard to verify the property 

D n (-Q) = {-l) n D n {Q) 



(13) 



which will be used in the forthcoming analysis and more precisely when we start to investigate the 
tunneling effect. Explicitly, the states 4>(x,y) are 



<f>(x, y) = [l B V^n\2 n ) a exp (-^\ H n (Q) 
where the Hermite polynomials H n {Q) are 

H n (Q) = (-l)"exp(Q 2 )^exp(-Q 2 ). 
Now injecting (fT2|) in Q, one can end up with 



2vlh 2 



,2 



l, vlh 2 



i 2 



D n (Q) = (e®) D n (Q) 



which is leading to the form 



2v 2 F h 2 
I 2 



V 



D n (Q) = (E®) D n (Q). 



It is clear that the eigenvalues read as 



a 2 n 



(14) 



(15) 



(16) 



(17) 



(18) 



where the constant a = —£ — . According to the interpretation of matter and antimatter, the bottom 

l B 

valley of the band structure takes the value — n that correspond to negative energies (anti- matter). 
For this reason, we write the eigenvalues as 



E$ = sgn(n)aV \n\ 
where the eigenfunctions are given by 



n E 



(19) 



(20) 



To complete our analysis, we need to determine the second spinor component, which can be 
obtained from the "kinetic balance" relation. Doing so to obtain 



1 



<Pi 



Em 



eB 



iv F I ip x +p y A x I D\ 



Jk y y 



It leads to the solution 



(f! = -sgn(n)i J D| n |_ 1 (Q)e 



ik y y 



(21) 



(22) 
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Finally, combining all to get the normalized eigenspinors as 

*<P k ,(x,y) = ±( -^ D H-,(x + x Q )\ lkyy 
{n ' k ^ ' V2\ D ]nl (x + x Q ) J ^ ! 

where s = sgn(ra) and the convention sgn(O) = should be taken into account. The energy spectrum 
has zero-mode wavefunction, such as 



</$ (a- y ) = _L [ ) (24) 

>,fc^ >"> { D (x + x ) ' 



These results are concerning the first region, which is assimilated to a subsystem of Dirac fermions 
in constant magnetic field. Next we will see how the above results will change when we move to the 
second region, which is also Dirac fermions but with a mass term. 

2.2 Second region spectrum 

From the nature of the system under consideration, we write the Hamiltonian corresponding to region 
(II) in terms of matrix as 

HM=v F ( ° «'-i*v) + (* V (25) 

\TT x + iTT y J \ -t' j 

In the Landau gauge, we have the form 

ffP) =v F { * R P *~ IP * :^ X ). (26) 

V Px + ipy + i^fx -V J 

Note that, the energy gap behaves like a mass term. Certainly this will affect the above results and 
lead to interesting consequences in underlying the basics features of such system. 

We need to derive the energy spectrum and its eigenspinors. In doing so, let us fix c/>( n ) = 
as a spmor of to write 



r'2 



which implies two relations 



Wp 



iv F (iPx + Py + ^-xj <p' 2 = - t') <p{ (28) 

iPx +Py + ^-x) v'l = (^ (H) + Pa- (29) 



These can be treated as we have done before to get one equation for one component spinor. After 
injecting (f28l) in (|29l) . we obtain a differential equation of second order for ip' 2 - This is 



h(p' 2 = 

It solution gives the second spinor component as 



<p' 2 . (30) 



x + x\ 



<f/ 2 (x, y) =D H ( ^— ^ ) e™, m£Z (31) 



where q y is a wave vector along y-direction in the second region and x' = q y l\- From last equation, 
it is easy to obtain the energy spectrum 

= sgrKWlyVH +t' 2 . (32) 

We notice that the term t' makes difference with respect to spectrum of region (I). It is convenient 
for our task to write energy as 



„/ r n , am 

ff( n ) -t'= —I 1 . 33 

TP. II) K I 



E<P) + 1> 

This will be used to show how the tunneling effect behaves in terms of the energy conservation and 
the parameter t' . 

Returning to the "kinetic balance" relations (j28[) and (I29p . one can easily derive the first spinor 
component (p±. It is 

-ia.\/\m\ _ / x + x' n 



M x >y) = i^^Ah> etw - (34) 



After normalization the eigenspinors read as 



^ Jx,y) = -1= [ -^M-l^ + ^o) ) e i W . (35) 



The normalization constants are given by 



, E( u ) +s't' Ei 11 ) -s't' 



(36) 



V EM ' V 
where we have set s' = sgn(m). They verify the useful relation 

a 2 m + b 2 m = 2. (37) 

We have also here a zero-mode wavefunction 

lU) V2\ D (x + x' ) J { ' 

Clearly, in the interface between two regions, there is conservation of the tangent components of the 
wave vector, i.e. k y = q y , which is equivalent to x$ = x' . This conclude the investigation of the energy 
spectrum for the present system in a constant magnetic field. Before using these results to solve some 
issues, it is interesting to start by underling properties of the eigenvalues. 



3 Selection rules 

Since we are considering two different regions characterized by two quantum numbers, one may inves- 
tigate their behaviors in terms of the energy system and the parameter if . For this, we use the energy 
conservation to establish a relation between them, which will be employed in the next in dealing with 
different issues. In particular, it will used to discuss the tunneling effect for our system through the 
reflexion and transmission coefficients. 
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Recall that, from the above analysis, we ended up with two energy spectra: and E^ l \ On 
the other hand, due to the system nature one should have an energy conservation, such as 



E 2 . 



(39) 



After replacing the energies by their expressions, it is easy to observe that the allowed energy values 
should verify the relation 

(40) 



t' 2 



\n\ — \m\ 

where the constraint \m\ < \n\ must be fulfilled. According to (|40p . one can realize n and m in terms 
of a pair of quantum numbers (p, q), such as 



±kp, 



m 



±k(p — q) 



(41) 



where | is an irreducible fraction with p > q and k is an integer value. This realization can be written 
in compact form as (n, m) 
values 



±k(p,p — q), which is equivalent to say that this pair lies in the set of 



{±(p,p- q), ±2(p,p - q),±3(p,p - q),±4(p,p -<?),-••}• 



(42) 



This is very important because without such set one can not talk about tunneling effect in the present 
case. We will clarify this statement from next section and exactly when we start to calculate different 
quantities in order to check the probability condition. 

The relation (|4(J[) is important in sense that its has some interesting consequences and can be used 
to give different interpretations of the present system. In fact, by inspecting some limiting cases, one 
can show 



• n = 

- E 2 



E 2 



m 



0. 



oo 



n 



On 



s, b % 



1. 



m 



0. 



In the next, we will see how these cases affect the tunneling effect through the evaluation of the 
reflexion and transmission coefficients. 

On the other hand, one can plot different quantum numbers to underline their behaviors. For this, 
we present three figures (1,2,3) illustrating some cases where the bold dots are exactly the allowed 
values for each numbers. They are 




Image001.eps 



7 



Figure 1: Variation of m in terms of re, it is showing that for n = 5 there are 6 possible values for m, 
i.e. m = 0, 1, • • • 5. 




Figure 2: Variation of the ratio of energy conservation and t', |, in terms of n. For n = 0, we have 
| = 1, which is in agreement with the first item cited above. 




Figure 3: Variation of the ratio of energy conservation and t' , |, in terms of m. For m = 0, we have 
identical behavior as for n = 0. 

4 Barrier in magnetic field 

To treat a concrete example of the present system, we consider the barrier in magnetic field. This 
latter can be seen as superposition of two regions separated by interface. We study the tunneling 
effect by evaluating the reflexion and transmission coefficients at interface, which corresponds to the 
point zero. These will be used to show that the probability condition is exactly one and emphasis 
what makes difference with respect to the same system studied in zero field [11J. To do this task, we 
distinguish two cases: propagation with positive and negative incidences. In both cases we deal with 
propagation from region (I) to (II) and vice verse. 

4.1 Propagation with positive incidence 

We start by defining what does mean propagation with positive incidence in our case. In fact, with 
this we will be able to do our job and therefore use the continuity equation at interface to determine 
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different quantities. Indeed, we split the obtained eigenspinors for two regions in positive and negative 
directions of the variable x, which is equivalent to write in region (I) 



,(i) = J_ -siD W _i(x + a? ) \ e%kyy Ai) = J_ siD\ n \-i(-x ~ x o) \ 
V2\ D H (x + x ) J V2 \ D H (-x-x ) J 

and in (II) we can do the same, such as 

^(11) _ 1 | -Omi D \m\-l( x + x o) \ e ik y y All) _ 1 ( - a miD\ m \_ x (-X - Xq) 
V%\ b m D\ m \(x + x ) J y/2 \ b m D\ m \(-x - x ) 



e ik y y ( 44 ) 



where we used the labels (±) to show which direction is concerned. We will see how these materials 
will serve as a tool to deal with different issues. For this, one has to consider two different cases, which 
will be the next task. 

Propagation from region (I) to region (II): This means that the waves are propagating from 
left to right and due to the continuity of the system, one can establish two equations. Indeed, at the 
interface x = and for all y, we can write 

^+r+ m ^=t+ r Jl l) (45) 

where r+ m and t+ m are reflexion and transmission coefficients, respectively, (+) indicates sign of 
propagation. From (03]) and (|4"1"|) . we end up with 

( -siD lnhl (x ) \ +r+ ( -siD[ n [_i(-a;o) \ =t + ( - a miD\ m \-i( x o) \ ^ 

nil x 0) / \ Um 

Using the property (fl~3j) . it is not hard to obtain 

r tm{- l )^ sD \n\-l{xo) + tn m a m D\m\-l( x o) = sD\ n \_x(x ) (47) 

-^n m (-l) |n| ^|n|(^o) + t+ m b m D\ m \(xo) = D H (x ). (48) 

They can be solved to explicitly get the coefficients in terms of some constants, which are magnetic 
field dependent. They are 

r+ _ (-l) |n| [sb m A nm (x ) - a m B nm (xo)] 
sb m A. nm [xo) -\- o, m B nm (^XQ) 

t+ m = 2sC n( x o) (5Q) 

sb m A nrn (xo^) + CL Tn B nm (^XQ^ 
where different involved quantities are given by 

Anm(xo) = D\ n \_ l (x )D\ m \(x ), B nm (x ) = D\ m \-l ( x o) D \n\ fco), 
C n (x ) = D\ nl _ 1 (x )D\ n \(x ). (51) 

These coefficients are just part of a set of quantities that should be completely determined. In fact, 
we still need other quantities to achieve our task, which can be obtained from next consideration. 

Propagation from region (II) to region (I): In similar way to the former case, we use the 
continuity at point zero to write 

€ ) +r+Jl I) =t+Jl\ (52) 



9 



Comparing this with ([45|) . we notice that there is an index interchange between different coefficients. 
By replacing, (f52|) becomes 



(54) 



( -a m «-D|m|-i(^o) \ +r+ ( -a m iD\ m \-i(- x o) \ _ f+ ( -isD\ n \_ x {x Q ) \ _ 
V KiD\ m \{x G ) J mn \ b m D\ m \(-xo) J mn \ D\ n \(x ) J 

These lead to the reflexion and transmission coefficients 

r+ (-l) im| [a m B nm (x ) - sb m A nm (x )] 

s bmA nm (xo) + a m i? nm (xo) 

^+ 2a m -F m (;co) (55) 

sb m A nm (xo) + o, m B nm {xo) 

where F m is given by 

F m {x ) = D\ m \_ 1 (x )D\ m \(x ). (56) 

This summarizes our analysis for propagation with positive incidence, which together will be used to 
discuss different issues and before doing so, we need to analyze negative incidence. 

4.2 Propagation with negative incidence 

One may ask about propagation with negative incidence. To reply this inquiry, we use the same 
analysis as before but one should take into account the negative sign of variable. For the neediness, 
we write the corresponding eigenspinors as 

$ = 2-1 - SW \n\-l(-* - *0) \ e - lkyy (I) = 1 ( -*iD H -l(* + *0) \ e - lky y (57) 

+ V2\ D H (-x-x ) J V2\ D H (x + x ) J 

and for second region, we have 

JIT) 1 I -a m iD\ m \_x(-x - x ) \ o _ ikyy jii) = _L_ ( -a m iD lmhl (x + x ) ^ ^_ ikyy 



V% \ b m D\ m \(-x -x ) J V2 \ b m D\ m \(x + x ) 



(58) 

By analogy to the former case, we distinguish the sense of propagation and determine the corresponding 
coefficients in terms of the magnetic field. 

Propagation from region (I) to region (II): Let us start from the continuity equation at point 
zero, such as 

^l ) +r- m ^ ) =t- m ^l I) (59) 



where label (— ) carried by the involved coefficients describes sense of propagation. Using (|57|) and 
to rewrite (1591 as 



( - siD \n\-l(-Xo) \ +r - ( - siD \n\-l(x ) \ = r ( - a ™iD\m\-l{-Xo) \ 
\ D \n\(-Zo) J D \n\{xo) ) b m D\ m \(-XQ) ) 

It can be solved to get 

(-l) |n| [sb m A nm (x ) - a m B nm (x )] 



sb m A nm {xo) -\- Q, m B nrn (xo) 



(61) 



t~ = (-l) |n|+H 2^n(x ) 

sb m A nm (xo) + Q>mB n m(xo) 
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Propagation from (II) to region (I): In similar way, we have at zero point 



*P + r m J U) = t-J? (63) 



showing the relation 



-a m iD| m |_i(-x ) \_ / -a m iD\m\-i(x ) \ = t - ( ~ siD M-i(-^o) \ ^ 
b m D\ m \{-xo) J mn \ b m D\ m \{x Q ) J mn \ D| n |(-x ) J 



The solutions are given by 



(-l)H [a m B nm (xo) - sb m A nm (x )] 
sb m A nm (x ) + a m B nm (x ) 



(65) 



t - = (-l)H+H 2 « m F m (x ) 

sb m A nm (xo) + a m B nm (xo) 

Finally we derived all materials needed to do our job. In fact, we use the obtained results to give 
different discussions and interpretations. 

4.3 Discussions 

We show the relevance of different coefficients derived so far. In fact, we will combine all of them 
to underline the basic features of the barrier in constant magnetic field. Certainly, from the above 
analysis one can notice that there are some relations between coefficients, which can be interpreted 
as a kind of symmetry. This will be helpful in sense that one can check the probability condition and 
derive other results. 

For the reason of do not repeating equations, let us choose a pair of index, such as i ^ j € {n, m}, 
to write them in compact forms. Indeed, from the above results, one can find 

To show the usefulness of these relations, let us define two coefficients in terms of the above quantities. 
The first one is 

p{xo) = ^-(xo)rS(xo) = — — — — 2 - (68) 

[sb m A nm (x ) + a m B nm (x )\ 



It is not hard to obtain 



rS(xo)rS(x ) = p(x )||^ = p(x )|^ = P (x ). (69) 

tij\ x 0) r ji\ x 0) 

In the same way we can write the second coefficient as 

/ \ .±, \.±, \ 4:sb m a m C n (xo)F m (x ) Asb m a m A nm (xo)B nm (x ) 

T{X0) = t^Xojt^Xo) = —— — — -—2 = — — — — — • (70) 

[sb m A nm (x ) + a m B nm (x )\ [sb m A nm (x ) + a m B nm (x )\ 
This implies that all quantities defined here can be related as 

C n (x )F rn (x ) = A nm (x )B nm (x ) (71) 
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which can easily check it from their expressions given before. Using (|68|) and (|7(jp . we show that they 
verify the condition 

p(x )+T(x ) = l. (72) 

This is among the interesting results derived so far. In fact, it tells us the transmission of barrier in 
magnetic field can not be greater than one, which is analogue to what obtained in zero field case [llj . 

At this level, one can inspect the relations (f68j) and (fTOj) to see how they behave with respect to 
three items cited in section 3. As well see soon these will change the nature of the present system. 
Indeed, taking into account such items, we straightforwardly end up with the results 

• n = => p(x Q ) = 0, t(xq) = 1. 

• fr — ► oo => p(x ) = 0, t(xq) = 1. 

• fr = 1 =>• p(xo) = 1, t(x ) = 0. 

These three cases have interesting interpretations in optics physics. More precisely, in the first and 
second points, the interface between two regions behaves like a non-reflective diopter, which means 
that everything is transmitted, i.e. total transmission. However, in last point, it looks like a mirror 
where the reflexion is total. 

To underline the basic properties of the transmission coefficient, we plot three figures in terms of 
the energy ratio E/t' where the conditions s = s' and E/t' ^ 1 are taking into account. These figures 
can be interpreted as follows. 

• Figure 4: One can see that each curve rises to an asymptote, which decreases when B nm /A nm 
increases. Clearly, for B nm /A nm = 1 the asymptote goes to 1, for B nm /A nm = 4 to 0.6 and 
finally for B nm /A nm = 8 to 0.4. 

• Figure 5: Each curve suddenly reached the value of 1 and then decreases to an asymptote, 
which decreases when B nm /A nm decreases. For B nm /A nm = 1 the asymptote goes to 1, for 
Bnm/A-nm = 1/4 to 0.6 and for B nm /A nm = 1/8 to 0.4. 

• Figure 6: Both curves tend towards the same asymptote for B nm /A nm = 1/2. The curve 
decreases to the asymptote 0.7 of the upper, but for B nm /A nm = 2 the curve rises to the 
asymptote of the lower. This behavior can be generalized to any case where B nm /A nm = I and 
B nm /A nm = 1/1, for all positive value I, which is in agreement with the statement cited in item 
2 that means that the transmission is total whenever ^ — > oo. 
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Figure 4: Variation of the transmission t(xq) in terms of the ratio E/t' for three cases: B nm = A nm , 
Bnm — 4j4 nm , B nm — 8j4 r 
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Figure 5: Variation of the transmission t(xq) in terms of the ratio E/t' but this case for B nm = A nm , 

Bnm — l/4j4ri}rt; B nm — l/8^4 nm . 
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Figure 6: Variation of the transmission t{xq) in terms of the ratio E/t' for two cases: B nm = 2A nm , 
Bnm — 1/2-Anm- 



5 Diode in magnetic field 

In last section we focused on the study of two different regions, with this it is natural to ask about a 
generalization to three regions where the second is characterized also with the energy gap t' and the 
first and third are identical. This is the case for instance of a diode in a constant magnetic field. To 
reply this inquiry, let us consider a system composed of a region indexed by the quantum number m 
of length 2w separating two others indexed by the same n. In fact, we will apply the same machinery 
as before to analyze the system behavior in such consideration. 

5.1 Reflexion and transmission coefficients 

The present situation is different from the former one, i.e. the barrier. This means that actually we 
have two interfaces fixed at the points wi = x$ — w and W2 = w + xo, which of course will make 
difference with respect of the barrier analysis. We use the above tool to write the continuity equation 
at each point and therefore derive the needed quantities for discussing tunneling effect in such case. 
In doing so, we consider also propagation with positive and negative incidences. 
In the first interface one can write 

-siD|„|_i(>i) \ +r+ f -si-D| n |-i(-^i) \ _ 
D\ n \(wi) J \ D|„|(-wi) J 
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= a( ~ a m iD \m\-l(wi) \ + p( - a m^| m |-l(-^l) \ ^ 

V b m D\ m \(wi) J \ b m D\ m \(-wi) J 

where r + , a and are three unknown parameters. The solution allows us to express two of them in 
terms of one for any point u>i = xq — w, such as 

a _ sbnAn^vm) + a m B nm {wi) + r + (-l)H [a m B nm {wi] - sb n A nm (wi)\ 

b m F m (wi) 

a _ a m Bnm(wi) - sb n A nm {wi) + r+(-l)l w l [a m B nm (wi) + sb n A nm (wi)] 

2{-l)\ m \a m b m F m {wi) 

These can be simplified to other equations by introducing two relevant relations in terms of w\. Indeed, 
using the former analysis one can show 

sb n A nm (wi) + a m B nm (wi) _ 1 
1d m b m F m 

(Wl) ttin{wi) 

a m B nm (wi) + sb n A nm (wi) = l ^\ n \+i r tm( w i) 

bmFm{w{) ttiniwi) 

With the help of these, we can rewrite a and [3 as 

a = - r + r im(wi) ^ 

tmn{wi) tmn{wi) 



/?=(-!) 



|n| + |m| 



(79) 



tmn(wi) tmn(wi) . 

As we will see, they will be employed to set up another couple of equation and therefore explicitly 
determine the reflexion and transmission coefficients for the diode case. 

To analyze the second interface, we start by considering the corresponding continuity equation at 
the point w 2 = xq + w. This is 

^ / a m iD\ m \_ x (w 2 ) \ + ^/_ 1 j| m | ( -a m iD\m\-i(w 2 ) \ = 
\ KiD\ m \(w 2 ) J \ b m D\ m \{w 2 ) J 



siD\ nl _ 1 (w 2 ) 

D\ n \(lV 2 ) 



(80) 



Note that, in obtaining (|80p in such form we made use of the property (|13p . The coefficients solutions 
are given by 

,+ sb m A nm (w 2 ) + a m B nm (w 2 ) t + 

OL — I — — — -~t — (Si) 

2am,bmFm{w 2 ) tmn{W 2 ) 



^ _ ^ a mB nr n(w 2 ) sb m A nm (w 2 ) _ ^ ^\n\ + \m\ r nm 



(w 2 ) 



2a m b m F m {w 2 ) ' ttin{w2) 

These lead to obtain 



(82) 



r + = r tm(wi) -r+ m (w 2 ) ^ 
1 - rt m {wi)rt m {w 2 ) 



From the above relations, one can find 



tmn{w\) 



1 ~ r nm( w l) r £m( w l) 

1 - rt m {wi)rtm{w 2 ) 



(84) 
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It convenient to write t + in terms of p and r as 



1 - p(wi) 



1 - r„ m (wi)r™(w2) 



*mn(wi) 



r(wi) 



1 - r„ m (roi)r„ m (w 2 ) 



(85) 



We close this part by nothing that the obtained results are only valid for the positive incidence. On 
the other hand, we will return to discuss the choice of fixing the points W\ and w 2 in terms of w. 

As far as the negative incidence is concerned, we can use a similar analysis to derive the corre- 
sponding coefficients. Indeed, from above we obtain 



r nm(»l)- r nm(« , 2) _ r+ m ( Wl ) - r+ m (w 2 ) 



1 - r nm (wi)r nm (w 2 ) 1 - r£ m (wi)rAm(w2) 



as well as 



tnmiwi) 



ThVl. 



1 - r n m(wi)r n m(w2) 



(W2 



t; 



(wi) 



t(wi) 



1 - rtm{wi)rnm{w2) 



(86) 



(87) 



These results will used to deal with the same issues as we have done for the barrier in magnetic field 
and also discuss other points. 



5.2 Collecting results 

After deriving the needed tools, we start by discussing the usefulness of them. For example, we can 
check the probability to see how the present system behaves at some critical points. In doing so, let 
us make an appropriate definition of the reflexion and transmission coefficients. It is convenient for 
our task to write 



R 



Explicitly, they take the forms 



R 



p(wi) + p(w 2 ) ~ 2r+ m (wi)r+ m (w 2 ) 
1 + p{wi)p{w 2 ) - 2rtm{wi)rt m {w2) 

1 - p(wi) - p(w 2 ) + p(w 1 )p(w 2 ) 
1 + p(wi)p(w 2 ) - 2rtm{wi)rtm{w 2 ) 



(89) 
(90) 



These show that how R and T for diode can be linked to the barrier coefficients at w\ and w 2 . 
Combining all to end up with probability one 



T + R = l. 



(91) 



This is among the interesting conclusion reached in this section, which analogue to that has been 
obtained for diode in zero field [TT] . 

One can see how the obtained results for the diode case will change when we consider our system 
centered around the point xq. In such case, its is not hard to show 



r nm(wi) = r+ m ( W2 ). 



(92) 



After injecting this relation in the forms of different coefficient derived here, one can conclude that 
everything is transmitted, i.e. the transmission is total. 
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5.3 Limiting case 

Another interesting case one should discuss is the resonant tunneling diode made of graphene in the 
presence of a magnetic field. This will allow us to characterize the system behavior and underline their 
properties. This limiting situation of the device is one where the barriers are described by a scalar 
Lorentz potential of the form [11] 

V(x,y) = limghl-e(\x\-e)]a z = ga z S(x). (93) 

The connection with the true barrier is made by identifying g with 2'yt'wa where 7 a numerical 
constant of dimensions inverse of length and a is the carbon-carbon distance. This form of the 
potential is equivalent to a mass term and therefore to a gap in the spectrum. However, given the 
short range nature of the potential, its effect comes only in the boundary conditions imposed on the 
wavefunction at the potential position. Note in passing that, the problem of Dirac electrons in delta 
function potentials has been studied in the past [14J and for a recent review one can consult |llj . 

To study (f93"l) in terms of our language, we adopt the same method as has been applied to the 
diode case for zero field |llj . It is based on the boundary condition around the point x = of the 
eigenspinors. This is equivalent to write 

(94) 

\ ¥>2(U' J J \ ¥>2(U j J 

where the matrix M is given by 

(cosh(g) isinh(g) \ , , 

• v,rt • (95) 

— zsmn(g) cosn(g) / 

The points ± represent positive and negative infinitesimals, the constant is g = ^ a . For some 
reasons, we need the inverse of (I94p . which can be obtained by determining M . It is easy to get 

(96) 




M -l _ ( cosh(g) -isinh(g) 
\ 2sinh(g) cosh(g) 



which together they verify the unitary condition M l M = I2. Combining all to end up with the 
inverse of (|M|) . such as 

After getting these two equations, one can use the former analysis to show that how the obtained 
results for diode in magnetic field will be affected. For this, we consider also two cases: propagation 
with positive and negative incidences. 

Propagation with positive incidence: Returning to our tools, we can map (|94p in terms of 
the cylindrical parabolic functions as 




-s^M-iO^o) ^(-I^'^h-iOeo) 
D H (x ) (-l)HZ) w (x ) 




(98) 
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c+ ( -a»J3| n |_i(a;o) 



(99) 



and in similar way, we have 

( <Pi (0+) ^ 

V ^(o+) t 

where A + ,B + and C + are three parameters those should be determined. According to (|94|) . the last 
equation takes the form 

( ^ l(0+) \=m( ~ SiD \n\-l^) si(-l)\ n \D {nhl (x ) 

{ J V D n( x o) (-ipD H (x ) ){b- 

This can be solved to obtain A + and B + as 



(100) 



A + = sC + 



2scosh.(g)D\ n \(x )D\ n \_ l (xo) + smh{g) D?,(x ) + Df._ ^Xq) 



B+ 



? C + (-l 



2 D \n\( x o) D \n\-l(xo) 
sinh(§) [Df^xo) - Df n[ (x ) 



(101) 



(102) 



2D| n |(ar )-D|n|-i(^o) 

The usefulness of these is to define the reflexion and transmission coefficients for the present situation. 
Indeed, After calculation, we show 



B+ 



(-l)H sinh(<?) 


D? 


_i(*o) - 






2scosh.(g)D\ n \(xo)D\ n \^ 1 (x ) + 


sinh(g) 


'Df n] (x )+Dl hl (x )_ 



2sD\ n \(x )D\ n \_ 1 (x ) 



t+ _ £1 = 

A+ 2scosh(g)D lnl (x )D lnl _ l (x ) + sinh(<?) D^(x ) + D^_ x (x Q ) 

Propagation with negative incidence: In this case, we can write for + 

( Vi(0 + ) \( siD lnhl (x ) si(-lpD lnhl (x ) \ ( A- \ 
\ M0 + ) J V D \^) (-lpD H (x ) 

as well as for 0~ 



(103) 
(104) 



B- 



I pi(o-) 
V ^2(0-) 

Using the same technique as above, one can obtain 
B~ 



-siD\ n \_x(x ) 



(-l)H sinh(ff) 






2s cosh(g)D\ n \(x )D\ n \_ 1 (x ) - sinh(g) 





f 



c- 



2sD\ n \(x )D\ n \_ l (x ) 



A 2scosh(5)Z?| n |(x )-D| n |_i(xo) - sinh(<?) Df^(x ) + Df nl l (x ) 



I' 



(105) 
(106) 

(107) 
(108) 



Having all ingredients, let us check the probability condition. This can be achieved by defining 
two quantities in terms of the above coefficients, such as 



R = r r 



-sinh 2 (<7) 


^(^-■Dfnl-lC^) 


2 


4cOsh 2 (5)Z)2 i| ( 2 ; )D2 j| _ 1 ( Xo ) _ sinh 2(£) 


5iN + ^ M wf 



4Df„[(«o)^|»|-i(^) 



4cOsh 2 (5) J D 2 n|(xo ) jD 2 i| _ i(:ro) _ sinh 2 ( - ) \ D fJ Xo)+D 2 ni _ i{xo) 



2 ■ 



(109) 
(110) 
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After a straightforward calculation, one can find a probability one 



R + T = 



4£f n| (x Q )L>f nM (x ) - sinh 2 (3) D* n \(x ) - D^xo) 




(111) 



4co S h 2 (5) J Df n| (xo) J Df n| _ 1 (^o) - sinh 2 Q) Df (x ) + Df (x ) 



This also has been obtained by studying the scalar Lorentz potential in zero field |llj . 

Let us analyze some limits of the transmission coefficient. This can be reached by inspecting (|109[) 
and as well as (|11U|) . Clearly, there is a trivial solution, i.e. T = 1, which can easily be obtained by 
requiring that the involved parameter g is nothing but zero. Furthermore, other interesting discussions 
can be reported by looking some cases involving the parameter g. After calculation, we end up with 
the results 

• g <C vph => T — > 1. 



Of course, one can also try to plot different figures in order to give a full description of the transmission 
behavior with respect to the scalar Lorentz potential taken here. 

6 Conclusion 

The present paper is devoted to the study of the tunneling effect for Dirac fermions in the presence 
of a constant magnetic field B. More precisely, we started by splitting our system in two different 
regions where the second one is characterized by an energy gap t' . Moreover, we derived the energy 
spectrum and the corresponding eigenspinors for each region. They are used to deal with our task 
and in particular show how Dirac fermions behave in such system. From the energy conservation E, 
we established an interesting relation between two quantum numbers (n, m) involved in game and the 
system energy. This allowed us to control the behavior of the energy and therefore fix the allowed 
values of (n, m) where their behaviors are plotted in the figures (1, 2, 3). In analyzing their properties, 
we noticed that there are three limiting cases those have significant influences on the study of the 
tunneling effect. 

Subsequently, we concentrated on two examples of the present system. Indeed, by considering a 
barrier in magnetic field, we analyzed the propagation for positive and negative incidences. In fact, 
using the continuity equations we determined the corresponding reflexion and transmission coefficients. 
They are obtained as a function of the ratio §■ and more importantly, they verified the probability 
condition, i.e. their sum is equal one. To characterize transmission of the system, we plotted the 
figures (4,5,6) showing its behaviors for different values of the involved parameters. Using the three 
limiting cases of (n, m), we noticed that the system can be seen either as a diopter or mirror depending 
to which case is included. These actually are corresponding to a total reflexion or total transmission, 
respectively. 

We focused on another interesting case that is a diode in magnetic field. This is equivalent to 
consider a system composed of three regions where the second regions is completely different from the 
first and third, which they are identical. After getting different coefficients for positive and negative 
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incidences, we defined two new quantities to end up with their sum is equal one, i.e. probability one. 
Requiring some constraint on the system, we showed that it is possible to have a total transmission. 
Moreover, we considered the case where the barriers are described by a scalar Lorentz potential of the 
form given in (|93p . This allowed us to study the resonant tunneling of diode make in graphene where 
different coefficients are obtained and their probability condition is verified. Moreover, we showed that 
there is a possibility to have a total transmission in such case. 

Finally, our findings so far will not remain at this stage. In fact, still some interesting questions one 
should answer to get a full descriptions of different issues. These concern for instance the Fabry-Perot 
solid etalon in the usual optics [15] and the related finesse factor. Since we have transmission, one can 
also ask about the conductivity to discuss the anomalous quantum Hall effect for the present system. 
Hopefully to come back to these questions in a subsequent work. 
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